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9.2 Frequency Response Characteristics 

The frequency transfer functions are defined for sinusoidal inputs having all possible 
frequencies uj £ [0, +oo). They are obtained from (9.1) by simply setting s = ju, 
that is 


Open—loop: G(ju)H(joj), Closed—loop: 


G(M 


1 + G(juj)H{juj) 


= M(ju) (9.1) 


Typical diagrams for the magnitude and phase of the open-loop frequency transfer 
function are presented in Figure 9.1. 
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Figure 9.1: Magnitude (a) and phase (b) of the open-loop transfer function 
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System Bandwidth: This represents the frequency range in which the magnitude 
of the closed-loop frequency transfer function drops no more than 3 dB (decibels) 
from its zero-frequency value. The system bandwidth can be obtained from the 
next equality, which indicates the attenuation of 3dB, as 


| AT (jujBw) | 


^=|M(0)| ujbw 


(9.2) 


Peak Resonance: This is obtained by finding the maximum of the function 
\M(ju)\ with respect to frequency uj. It is interesting to point out that the 
systems having large maximum overshoot have also large peak resonance. This 
is analytically justified for a second-order system in Problem 9.1. 

Resonant Frequency: This is the frequency at which the peak resonance occurs. 
It can be obtained from 


l \ M U .)\=0 


(jJf 



Figure 9.2: Magnitude of the closed-loop transfer function 
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9.3 Bode Diagrams 

Bode diagrams represent the frequency plots of the magnitude and phase of the 
open-loop frequency transfer function G{ju)H{ju). The magnitude is plotted in 
dB (decibels) on the log a; scale. We first study independently the magnitude and 
frequency plots of each of these elementary frequency transfer functions. Since the 
open-loop frequency transfer function G{ju)H{ju) is given in terms of products 
and ratios of elementary transfer functions, it is easy to see that the phase of 
G{ju)H{ju) is obtained by summing and subtracting phases of the elementary 
transfer functions. Also, by expressing the magnitude of the open-loop transfer 
function in decibels, the magnitude \G{ju)H{ju)\ dB is obtained by adding the 
magnitudes of the elementary frequency transfer functions. For example 


\G{ju)H{ju)\ dB = 201og 10 


K{ju + z\){ju + zf) 
{ju){ju + pf){ju + p%) 


201og 10 \K\+ 201og 10 | ju + zi\ + 201og 10 | ju + Z2 


+ 20 log 10 


ju 


+ 201og 10 


ju + P 2 


+ 20 log 10 


ju + P3 


and 


arg {G(ju)H {ju)} = arg {K} + arg {ju + z\} + arg {ju + zf} 

- arg {ju} - arg {ju + p 2 } - arg {ju + p 3 } 
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Constant Term: Since 

K c ib = 20 log!o K 


positive number K > 1 
negative number K < 1 


f 0°, K > 0 
\ —180°, K< 0 


(9.3) 


the magnitude and phase of this element are easily drawn and are presented in 
Figure 9.3. 



Figure 9.3: Magnitude and phase diagrams for a constant 


Pure Integrator: The transfer function of a pure integrator, given by 

G(jcu) = — (9.4) 

has the following magnitude and phase 

I G(ju)\ dB = 201og 10 - = —201og 10 to, arg G(ju) = -90° (9.5) 

UJ 
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It can be observed that the phase for a pure integrator is constant, whereas the 
magnitude is represented by a straight line intersecting the frequency axis at uj — 1 
and having the slope of — 20dB/decade. Both diagrams are represented in Figure 
9.4. Thus, a pure integrator introduces a phase shift of — 90° and a gain attenuation 
of — 20dB/decade. 
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Figure 9.4: Magnitude and phase diagrams for a pure integrator 
Pure Differentiator: The transfer function of a pure differentiator is given by 

G{ju) = ju (9.6) 

Its magnitude and phase are easily obtained as 

\G(ju)\ dB = 201ogio arg G(ju) = 90° (9.7) 







386 


FREQUENCY DOMAIN CONTROLLER DESIGN 


The corresponding frequency diagrams are presented in Figure 9.5. It can be 
concluded that a pure differentiator introduces a positive phase shift of 90° and 
an amplification of 20dB/decade. 

Real Pole: The transfer function of a real pole, given by 


G{jw) 


P 1 

p + juj 1+j- 


(9.8) 


has the following magnitude and phase 



Figure 9.5: Magnitude and phase diagrams for a pure differentiator 

r / \T /2 / \ 

\ G 0 w )\dB = -201°gio 1+ (“) > argG(jo;) =-tail -1 f (9.9) 
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The phase diagram for a real pole can be plotted directly from (9.10). It can 
be seen that for large values of uj, uj p , the phase contribution is —90°. For uj 
small, uj <C p, the phase is close to zero, and for uj = p the phase contribution is 
—45°. This information is sufficient to sketch arg G(juj) as given in Figure 9.6. 

For the magnitude, we see from (9.10) that for small uj the magnitude is 
very close to zero. For large values of uj we can neglect 1 compared to uj/p 
so that we have a similar result as for a pure integrator, i.e. we obtain an 
attenuation of 20dB/decade. For small and large frequencies we have straight- 
line approximations. These straight lines intersect at uj = p, which is also known 
as a corner frequency. The actual magnitude curve is below the straightline 
approximations. It has the biggest deviation from the asymptotes at the corner 
frequency (see Figure 9.6). 




Figure 9.6: Magnitude and phase diagrams for a real pole 
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Real Zero: The transfer function of an element representing a real zero is 

G(juj) = -(z + ju) = 1 + .7 (“) (9-10) 

Its magnitude and phase are 

r 21 1 / 2 

1^0’^)!^ = 201 °gio 1+ (“) ’ ^gG(juj) = tail -1 (^) (9.11) 

For small frequencies an asymptote for the magnitude is equal to zero and for large 
frequencies the magnitude asymptote has a slope of 20dB/decade and intersects the 
real axis at uj — z (the corner frequency). The phase diagram for small frequencies 
also has an asymptote equal to zero and for large frequencies an asymptote of 90°. 



Figure 9.7: Magnitude and phase diagrams for a real zero 
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Complex Conjugate Poles: The transfer function is given by 


G(jlo) 


w 


n 


(jujf + 2 (un(ju) + fl - 4) + .?2C— 

V W nj J 

The magnitude and phase of this second-order system are given by 

.2 \ 21 1/2 


\ G (ju)\ d B = —201ogi 0 


2<> 


+ 11-^ 


LJ. 


n 


LJ, 


n . 


arg G(juj) = — tan 


-l 


( 2 C,LO n 


v 4 - 


LJ 


LJ< 


For large values of u the corresponding approximations of (9.14) are 

. 2 ' 


\G(M\dB « —20log 10 


LJ 


LJ, 


—40 log 10 


LJ 


n. 


LJ\ 


n 


cirg{G(jLj)} — tan 1 ( ^ u ^ » — tan 1 (0 ) = —180 c 


—LJ 


(9.12) 


(9.13) 


At low frequencies the approximations can be obtained directly from (9.13), that is 

G(ju))tt- f = l =4> \G(jLj)\ dB = 0, arg {G(ju)} = 0° 

The asymptotes for small and large frequencies are, respectively, zero and 
—40dB/decade (with the corner frequency at lj = lj u ) for the magnitude, and 
zero and —180° for the phase. At the corner frequency u n the phase is equal to 
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—90°. Note that the actual plot in the neighborhood of the corner frequency depends 
on the values of the damping ratio £. Several curves are shown for 0.1 <c<i. 
It can be seen from Figure 9.8 that the smaller the higher peak. 



Figure 9.8: Magnitude and phase diagrams for complex conjugate poles 


Complex Conjugate Zeros: The complex conjugate zero is given by 

/ \ / • \ 2 / \ 2 / \ 


— 1 + 2(j ( ) + ( ) — 1 — ( ) + j%C ( 


(9.14) 


so that the corresponding Bode diagrams will be the mirror images of the Bode 
diagrams obtained for the complex conjugate poles represented by (9.13). Here, the 


asymptotes for small frequencies are equal to zero for both the magnitude and phase 
plots; for high frequencies the magnitude asymptote has a slope of 40dB/decade 
and starts at the corner frequency of u = u n , and the phase plot asymptote is 180°. 






FREQUENCY DOMAIN CONTROLLER DESIGN 


391 


9.3.1 Phase and Gain Stability Margins from Bode Diagrams 

Bearing in mind the definition of the phase and gain stability margins given in (4.54) 
and (4.55), and the corresponding phase and gain crossover frequencies defined in 
(4.56) and (4.57), it is easy to conclude that these margins can be found from Bode 
diagrams as indicated in Figure 9.9. 



Figure 9.9: Gain and phase margins and Bode diagrams 
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Example 9.1: Consider theopen-loop frequency transfer function 

{ju + 1 ) 


G(juj)H(juj) 


+ 2) (juj) + 2 (jto) + 2 


Using MATLAB 

num=[l 1];dl=[1 0];d2=[l 2];d3=[l 2 2]; 

denl = conv(dl, d2) ; den=conv(deni , d3); 

bode(num,den);[Gm,Pm,wcp,wcg]=margin(num,den); 


The phase and gain stability margins and the crossover frequencies are 
Gm = 8.9443dB, Pm = 82.2462°, u cp = 1.7989rad/s, u cg = 0.2558rad/s 

o 
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9.3.2 Steady State Errors and Bode Diagrams 

Steady state errors can be indirectly determined from Bode diagrams by reading 
the values for constants K p ,K v ,K a from them. Knowing these constants, the 
corresponding errors are easily found by using formulas (6.30), (6.32), and (6.34). 
The steady state errors and corresponding constants K p . K n .K a are first of all 
determined by the system type, which represents the multiplicity of the pole at the 
origin of the open-loop feedback transfer function, in general, represented by 


G{ju)H{ju)) 


K{ju + z\){ju + z 2 ) • • • 

(ju) r (ju + Pl)(ju + P2) 


This can be rewritten as 


where 


G(ju)H(ju) 


KziZ2 -(i + f)(i + f).- 
PlP2---W'(l + f)(l+f) ••• 


Kb {l + f)(l + f) ••• 

W'( 1 + f)( 1+ f) -• 


K b 


Kz\Z2 • • • 
PlP2 * * * 


(9.15) 


(9.16) 


(9.17) 


is known as Bode’s gain, and r is the type of feedback control system 



394 


FREQUENCY DOMAIN CONTROLLER DESIGN 


For control systems of type r — 0, the position constant according to formula 
(6.31) is obtained from (9.17) as 


Kb {l + f)(l + f) ••• 

(i-)°(i+ £)(! + £) - I* 


(9.18) 


It follows from (9.17)—(9.19) that the corresponding magnitude Bode diagram of 
type zero control systems for small values of u is flat (has a slope of 0 dB) and the 
value of 20 log Kb = 20 log K p . This is graphically represented in Figure 9.12. 


\G(j(n)H(j(o)\ 



Figure 9.12: Magnitude Bode diagram of type zero control systems at small frequencies 

For control systems of type r = 1, the open-loop frequency transfer function is 
approximated at low frequencies by 


i0(l + f)(l + f)--- 

w 1 ( 1 +^)( 1 +^)- 


K b 

( i ^) 1 


(9.19) 


It follows that the corresponding magnitude Bode diagram of type one control 
systems for small values of w has a slope of — 20dB/decade and the values of 


20 log 


Kb 

jw 


20log \ Kb\ — 20log | uj 


(9.20) 
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From (9.20) and (6.33) it is easy to conclude that for type one control systems the 
velocity constant is K v = Kb- Using this fact and the frequency plot of (9.21), we 
conclude that K v is equal to the frequency uo* at which the line (9.21) intersects 
the frequency axis, that is 

0 = 20log \Kb\ — 20log |u/| =>- Kb = a;* = K v (9.21) 

This is graphically represented in Figure 9.13. 


\G(j(o)H(j<o)\ dB \G(ja»H(jw)\ dB 



Figure 9.13: Magnitude Bode diagram of type one control systems at small frequencies 

Note that if K v = u* > 1, the corresponding frequency u* is obtained at the point 
where the extended initial curve, which has a slope of — 20dB/decade, intersects 
the frequency axis (see Figure 9.13b). 

Similarly, for type two control systems, r = 2, we have at low frequencies 

W 2 (i + ff)(i + ff) ••• ~ ( M 2 


( 9 . 22 ) 
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which indicates an initial slope of —40dB/decade and a frequency approximation 
of 


20 log 


K b 

(jw) 2 


20\og \Kb\ —20log (jj 2 


20log \Kb\ — 40 log a; 


(9.23) 


From (9.23) and (6.35) it is easy to conclude that for type two control systems the 
acceleration constant is K a = Kb- From the frequency plot of the straight line 
(9.24), it follows that Kb = (a/*) 2 , where u** represents the intersection of the 
initial magnitude Bode plot with the frequency axis as represented in Figure 9.14. 



Figure 9.14: Magnitude Bode diagram of type two control systems at small frequencies 


It can be seen from Figures 9.12-9.14 that by increasing the values for the 
magnitude Bode diagrams at low frequencies (i.e. by increasing Kb), the constants 
K P ,K V , and K a are increased. According to the formulas for steady state errors, 
given in (6.30), (6.32), and (6.34) as 


'SS 


step 


1 + K P ' 


-SS 


ramp 


K r 


'SSparabalic 


K a 
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we conclude that in this case the steady state errors are decreased. Thus, the bigger 
Kb, the smaller the steady state errors. 

Example 9.3: Consider Bode diagrams obtained in Examples 9.1 and 9.2. The 
Bode diagram in Figure 9.10 has an initial slope of —20 dB/decade which intersects 
the frequency axis at roughly cu* = 0.2 rad/s. Thus, we have for the Bode diagram 
in Figure 9.10 


K p = oo, K v « 0.2, K a = 0 


Using the exact formula for K v , given by (6.33), we get 

(s+l) 


K, = lim < s- 


= 0.25 


>6 ^ s(s + 2)(s 2 + 2s + 2) 

In Figure 9.11 the initial slope is OdB, and hence we have from this diagram 


20 log K p ~ 15 =>- Kp ~ 5.62, K v = 0, K a = 0 


Using the exact formula for K p as given by (6.31) produces 


K t) = lim 

1 s -^0 


10(s + 4) 


(s + l)(s + 2)(s + 3) 


6.67 


Note that the accurate results about steady state error constants are obtained easily 
by using the corresponding formulas; hence the Bode diagrams are used only for 
quick and rough estimates of these constants. 


o 
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9.4 Compensator Design Using Bode Diagrams 

Controller design techniques in the frequency domain will he governed by the 
following facts: 

(a) Steady state errors are improved by increasing Bode’s gain Kb- 

(b) System stability is improved by increasing phase and gain margins. 

(c) Overshoot is reduced by increasing the phase stability margin. 

(d) Rise time is reduced by increasing the system’s bandwidth. 

The first two items, (a) and (b), have been already clarified. In order to justify 
item (c), we consider the open-loop transfer function of a second-order system 


n 


(jw)(jw + 2(co n ) 

whose gain crossover frequency can be easily found from 

\G(ju cg )H(ju cg )\ = , ^ . = 1 


(9.24) 


UJ^UJ 2 Tipcul 


(9.25) 


leading to 


U C g = UJ n J \ A+ 4 C 2 - 2 C 2 


The phase of (9.25) at the gain crossover frequency is 

arg {G(ju C g)H(ju cg )} = -90° - tan 


-l 


u, 


c 9 


2 C,uJn 


so that the corresponding phase margin becomes 

2 C 


Pm = tan" 


V 7 ! + 4C 2 - 2C 2 


Pm(C) 


(9.26) 


(9.27) 


(9.28) 
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Plotting the function it can be shown that it is a monotonically increasing 

function with respect to £; we therefore conclude that the higher phase margin, the 
larger the damping ratio, which implies the smaller the overshoot. 

Item (d) cannot be analytically justified since we do not have an analytical 
expression for the response rise time. However, it is very well known from 
undergraduate courses on linear systems and signals that rapidly changing signals 
have a wide bandwidth. Thus, systems that are able to accommodate fast signals 
must have a wide bandwidth. 


9.4.1 Phase-Lag Controller 

The transfer function of a phase-lag controller is given by 

f)f\z[ + ;juj 1+i^ 




- 


z\ J Pi + ju) 1 + j 


iO 5 

Pi 


Zl > Pi 



(9.29) 


Figure 9.15: Magnitude approximation and exact phase of a phase-lag controller 
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Due to attenuation of the phase-lag controller at high frequencies, the frequency 
bandwidth of the compensated system (controller and system in series) is reduced. 
Thus, the phase-lag controllers are used in order to decrease the system bandwidth 
(to slow down the system response). In addition, they can be used to improve the 
stability margins (phase and gain) while keeping the steady state errors constant. 

Expressions for uj max and 4> max of a phase-lag controller will be derived in the 
next subsection in the context of the study of a phase-lead controller. 


9.4.2 Phase-Lead Controller 

The transfer function of a phase-lead controller is 


Glead(j(d) 


Pi \ z -2 + ju) _ 1 +Jjr 2 
Z2jP2+ju 


P2 > Z2 


(9.30) 



Figure 9.16: Magnitude approximation and exact phase of a phase-lead controller 
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Due to phase-lead controller (compensator) amplification at higher frequencies, it 
increases the bandwidth of the compensated system. The phase-lead controllers 
are used to improve the gain and phase stability margins and to increase the system 
bandwidth (decrease the system response rise time). 

It follows from (9.31) that the phase of a phase-lead controller is given by 


-l { ^ \ / U) 


arg {Gi eci d(j^)} = tan 1 I — I - tan 

Z2j \P2 


so that 


d_ 

dw 


arg — 0 to 

max = \J Z‘2P‘2 


Assume that 


P 2 = az 2 , a > 1 


^ max — 


P2 


Substituting u max in (9.32) implies 


tail 4>max — 


a — 1 
2 y/a 


The parameter a can be found from (9.35) in terms of 4>max as 

1 + sin cf), 


a = 


mi a x 


1 — sin (f>, 


max 


(9.31) 

(9.32) 

(9.33) 

(9.34) 

(9.35) 


Note that the same formulas for uj max , (9.33), and the parameter a, (9.36), hold 
for a phase-lag controller with pi, z\ replacing p 2 ,Z 2 and with p\ = az\, a < 1. 
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9.4.3 Phase-Lag-Lead Controller 

The phase-lag-lead controller has the features of both phase-lag and phase-lead 
controllers and can be used to improve both the transient response and steady state 
errors. The frequency transfer function of the phase-lag-lead controller is given by 


ri ( ■ \ {ju + Zl)(jLJ + Z 2 ) ZlZ2\ 1+ ^z 1 )[ 1 ^^z 2 ) 

G cljUJ) = 7T-—-r =-7-r—7-r- 

(JU+P1)(JU+P2) P1P2 M + j!£.\ 


(i + j1)(i+^) 

( 1 + -?'b)( 1 +j k) 


Z\Z 2 = piP2 1 P2> Z2> ZI> pi 


The Bode diagrams of this controller are shown in Figure 9.17. 


(9.36) 



Figure 9.17: Bode diagrams of a phase-lag-lead controller 
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9.4.4 Compensator Design with Phase-Lead Controller 

The following algorithm can be used to design a controller (compensator) with a 
phase-lead network. 

Algorithm 9.1: 


1. Determine the value of the Bode gain Kb given by (9.18) as 

Kz\Z2 • • • 

Kb =- 

P1P2 • • • 

such that the steady state error requirement is satisfied. 

2. Find the phase and gain margins of the original system with Kb determined 
in step 1. 

3. Find the phase difference, A</>, between the actual and desired phase margins 
and take (f) max to be 5°-10° greater than this difference. Only in rare cases 
should this be greater than 10°. This is due to the fact that we have to give an 
estimate of a new gain crossover frequency, which can not be determined very 
accurately (see step 5). 

4. Calculate the value for parameter a from formula (9.36), i.e. by using 


1 + sin cfimax 
1 sill <f>max 


> 1 


5. Estimate a value for a compensator’s pole such that cj max is roughly located at 
the new gain crossover frequency, Lo max ~ a; C gnew As a rule of thumb, add the 
gain of AG = 201og(n)[dB] at high frequencies to the uncompensated system 
and estimate the intersection of the magnitude diagram with the frequency axis, 
say uj i. The new gain crossover frequency is somewhere in between the old 
uj cg and uj [. Some authors (Kuo, 1991) suggest fixing the new gain crossover 
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frequency at the point where the magnitude Bode diagram has the value of 
—0.5AG[dB]. Using the value for parameter a obtained in step 4 find the value 
for the compensator pole from (9.34) as — p c = —u max \fa and the value for 
compensator’s zero as —z c = —p c /a. Note that one can also guess a value for 
p c and then evaluate z c and u max . The phase-lead compensator now can be 
represented by 



as + p c 
s + Pc 


6. Draw the Bode diagram of the given system with controller and check the values 
for the gain and phase margins. If they are satisfactory, the controller design is 
done, otherwise repeat steps 1-5. 

Example 9.4: Consider the following open-loop frequency transfer function 


G(ju)H(ju) 


_ K{ju + 6) _ 

(jcu + + 2 )(jto + 3) 


Step 1. Let the design requirements be set such that the steady state error due to a 
unit step is less than 2% and the phase margin is at least 45°. Since 


1 

1 + K p 


1 

1 + K b ' 


K b 


K x 6 
1x2x3 


= K 


we conclude that K > 50 will satisfy the steady state error requirement of being 
less than 2%. We know from the root locus technique that high static gains can 
damage system stability, and so for the rest of this design problem we take K = 50. 
Step 2. We draw Bode diagrams of the uncompensated system with the Bode gain 
obtained in step 1 and determine the phase and gain margins and the crossover 
frequencies. This can be done via MATLAB 
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[den] = [16 11 6] ; 

[num] = [50 3 0 0]; 

[Gm,Pm,wcp,wcg]=margin(num,den); 
bode(num,den) 

The corresponding Bode diagrams are presented in Figure 9.18a. The phase and 
gain margins are obtained as Gm = oo, Pm = 5.59° and the crossover frequencies 
are uj cg = 7.5423 rad/s, u cp = oo. 

Step 3. Since the desired phase is well above the actual one, the phase-lead controller 
must make up for 45° — 5.59° = 39.41°. We add 10°, for the reason explained in 
step 3 of Algorithm 9.1, so that 4> max = 49.41°. 



Figure 9.18: Bode diagrams for the original system 
(a) and compensated system (b) of Example 9.4 
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The above operations can be achieved by using MATLAB as follows 

% estimate Phimax with Pmd = desired phase margin; 
Pmd=input('enter desired value for phase margin'); 
Phimax=Pmd-Pm+l0; 

% converts Phimax into radians; 

Phirad=(Phimax/180)*pi; 

Step 4. Here we evaluate the parameter a according to the formula (9.36) and get 
a = 7.3144. This can be done in MATLAB by 

a= (1 + sin(Phirad) )/(1-sin(Phirad) ) ; 

Step 5. In order to obtain an estimate for the new gain crossover frequency 
we first find the controller amplification at high frequencies, which is equal to 
20 log (a) = 17.2836 dB = A G c ib- The magnitude Bode diagram increased by 
A G ( ib at high frequencies intersects the frequency axis at u ma x ~ 10.5 rad/s. 
We guess (estimate) the value for p c as p c = 25, which is roughly equal to 
u ma x\/u- By using p c = 25 and forming the corresponding compensator, we get 
for the compensated system Pmc = 48.2891° at u cgnew = 13.8519 rad/s, which is 
satisfactory. This step can be performed by MATLAB as follows. 

% Find amplification at high frequencies, DG; 
DG=20*logl0(a); 

% estimate value for pole —pc from Step 5; 
pc=input('enter estimated value for pole pc'); 

% form compensator's numerator; 
nc=[a pc]; 

% form compensator's denominator; 
dc=[1 pc]; 

% find the compensated system transfer function; 
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numc=conv(num,nc) ; 
denc=conv(den,dc) ; 

[Gmc,Pmc,wcp,wcg]=margin(numc, dene) ; 
bode(numc,dene) 

The phase-lead compensator obtained is given by 

, . 7.3144s T 25 as P 2 
K) s + 25 s T p 2 

Step 6. The Bode diagrams of the compensated control system are presented in 
Figure 9.18b. Both requirements are satisfied, and therefore the controller design 
procedure is successfully completed. 

Note that num, den, numc, dene represent, respectively, the numerators and 
denominators of the open-loop transfer functions of the original and compensated 
systems. In order to find the corresponding closed-loop transfer functions, we use 
the MATLAB function cloop, that is 

[enum,eden]=cloop(num,den,-1); 

% —1 indicates a negative unit feedback 
[enume,cdenc]=cloop(numc, dene, -1); 

The closed-loop step responses are obtained by 

[y,x]=step(enum,eden); 

[yc,xc]=step(enume,cdenc); 

and are represented in Figure 9.19. It can be seen from this figure that both 
the maximum percent overshoot and the settling time are drastically reduced. In 
addition, the rise time of the compensated system is shortened since the phase-lead 
controller increases the frequency bandwidth of the system. 
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Figure 9.19: Step responses for the original (a) and compensated (b) systems 



Figure 9.20: Root locus of the original system 
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9.4.5 Compensator Design with Phase-Lag Controller 

Compensator design using phase-lag controllers is based on the compensator’s 
attenuation at high frequencies, which causes a shift of the gain crossover frequency 
to the lower frequency region where the phase margin is high. The phase-lag 
compensator can be designed by the following algorithm. 

Algorithm 9.2: 

1. Determine the value of the Bode gain Kb that satisfies the steady state error 
requirement. 

2. Find on the phase Bode plot the frequency which has the phase margin equal to 
the desired phase margin increased by 5° to 10°. This frequency represents the 
new gain crossover frequency, u cg new 

3. Read the required attenuation at the new gain crossover frequency, i.e. 
\AG(juj cgnew )\ dB , and find the parameter a from 

-20log = —20log (a) = \AG(jcj cgnew )\ dB 


which implies 


Note that 


a = 




cgnewJlfiB 


1 

I AG(juj cgn eie) | 


\AG(jLJcgn ew ) | 


\K\\j^cgnem T - j^cg new + 2:2 
| j^cgnew Pi j^cgnew + P 2 | * 
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4. Place the controller zero one decade to the left of the new gain crossover 
frequency, that is 


0}eg new 


Find the pole location from p c = az c = au cgnew / 10. The required compensator 
has the form 



as + p c 
5 + Pc 


5. Redraw the Bode diagram of the given system with the controller and check the 
values for the gain and phase margins. If they are satisfactory, the controller 
design is done, otherwise repeat steps 1-5. 


Example 9.5: Consider a control system represented by 

G(s \ = K 
K 1 s(s + 2 )(s + 30) 

Design a phase-lag compensator such that the following specifications are met: 
e SSr < 0.05, Pm > 45°. The minimum value for the static gain that produces 
the required steady state error is equal to K = 1200. The original system with this 
static gain has phase and gain margins given by Pm = 6.6449°, Gm = 4.0824 dB 
and crossover frequencies of uj cg = 6.1031 rad/s, uj cp = 7.746 rad/s. 

The new gain crossover frequency can be estimated as uj cgn ew = 1.4 rad/s since 
for that frequency the phase margin of the original system is approximately 50°. At 
ujegnew = 1.4 rad/s the required gain attenuation is obtained by MATLAB as 
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wcgnew=l.4; 
dl = 12 0 0; 
gl=abs(j *1.4) ; 
g2=abs(j *1.4 + 2); 
g3=abs(j *1.4 + 30); 
dG=dl/(gl*g2 *g3); 

which produces |AG(jl.4)| = 11.6906 and a = l/|AG(jl.4)| = 0.0855. The 
compensator’s pole and zero are obtained as —z c = —uJcgnew/ 10 = —0.14 and 
— p c = —aucgnew/ 10 = —0.0120 (see step 4 of Algorithm 9.2). The transfer 
function of the phase-lag compensator is 



0.0855s+ 0.0120 
s +0.0120 



10 -1 10° 10 1 io 2 

Frequency (rad/sec) 


Figure 9.21: Bode diagrams for the original system 
(a) and compensated system (b) of Example 9.5 
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The new phase and gain margins and the actual crossover frequencies are Pmc = 
47.03°, Gmc = 24.82 dB, uj C gnew = 1.405 rad/s, uicpnew = 7.477rad/s and 
so the design requirements are satisfied. The step responses of the original and 
compensated systems are presented in Figure 9.22. 



Figure 9.22: Step responses for the original system 
(a) and compensated system (b) of Example 9.5 

It can be seen from this figure that the overshoot is reduced from roughly 0.83 to 
0.3. In addition, it can be observed that the settling time is also reduced. Note that 
the phase-lag controller reduces the system bandwidth (uicgnew < u cg ) so that the 
rise time of the compensated system is increased. 


o 
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9.4.6 Compensator Design with Phase-Lag-Lead Controller 

Compensator design using a phase-lag-lead controller can be performed according 
to the algorithm given below, in which we first form a phase-lead compensator and 
then a phase-lag compensator. Finally, we connect them together in series. 

Algorithm 9.3: 

1. Set a value for the static gain Kb such that the steady state error requirement 
is satisfied. 

2. Draw Bode diagrams with Kb obtained in step 1 and find the corresponding 
phase and gain margins. 

3. Find the difference between the actual and desired phase margins, 
A <p = Pmd — Pm, and take p max to be a little bit greater than A <p. Cal¬ 
culate the parameter a 2 of a phase-lead controller by using formula (9.36), 
that is 

1 + sin (f>max 

a '2 = 1 - • , 

1 Sill <pmax 

4. Locate the new gain crossover frequency at the point where 

20 log \ G(j 

^ eg new )| = —10 log n 2 (9.37) 

5. Compute the values for the phase-lead compensator’s pole and zero from 

Pc2 ~ ^cgnew\ /^2 1 Zc'2 = Pc2/&2 (9.38) 

6. Select the phase-lag compensator’s zero and pole according to 


Z c l = 0.1z c2 , Pci = Zc\ /«2 


( 9 . 39 ) 
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7. Form the transfer function of the phase-lag-lead compensator as 

ri ( \ n / \ / \ S -\- Z c \ S + Z C 2 

G c (s) — Gi a g(s) X Gi ea d(s) — -;-X —■- 

5+ Pd S+p c 2 

8. Plot Bode diagrams of the compensated system and check whether the design 
specifications are met. If not, repeat some of the steps of the proposed algo¬ 
rithm—in most cases go back to steps 3 or 4. 

The phase-lead part of this compensator helps to increase the phase margin 
(increases the damping ratio, which reduces the maximum percent overshoot and 
settling time) and broaden the system’s bandwidth (reduces the rise time). The 
phase-lag part, on the other hand, helps to improve the steady state errors. 

Example 9.6: Consider a control system that has the open-loop transfer function 

C(s) = K ^ S - 10 ) 

K ’ (s 2 + 2s + 2)(s + 20) 

For this system we design a phase-lag-lead controller by following Algorithm 9.3 
such that the compensated system has a steady state error of less than 4% and 
a phase margin greater than 50°. In the first step, we choose a value for the 
static gain K that produces the desired steady state error. It is easy to check that 
K = 100 =>- e ss = 3.85%, and therefore in the following we stick with this value 
for the static gain. Bode diagrams of the original system with K = 100 are pre¬ 
sented in Figure 9.23. 
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Gm=lnf dB, (w= NaN) Pm=31.61 deg. (w=7.668) 



Frequency (rad/sec) 


Figure 9.23: Bode diagrams of the original system 

It can be seen from these diagrams—and with help of MATLAB determined accu¬ 
rately—that the phase and gain margins and the corresponding crossover frequen¬ 
cies are given by Pm = 31.61 °,Gm = oo, and uo cg = 7.668rad/s, uj cp = oo. 
According to step 3 of Algorithm 9.3, a controller has to introduce a phase lead 
of 18.39°. We take 4> m ax = 25° and find the required parameter a 2 = 2.4639. 
Taking u cgnew = 20 rad/s in step 4 and completing the design steps 5-8 we find 
that Pm = 39.94°, which is not satisfactory. We go back to step 3 and take 
dntax = 30° = 0.5236 rad, which implies a -2 = 3. 

Step 4 of Algorithm 9.3 can be executed efficiently by MATLAB by performing 
the following search. Since —10 log 3 = —10.9861 dB we search the magnitude 
diagram for the frequency where the attenuation is approximately equal to —11 dB. 
We start search at u = 20 rad/s since at that point, according to Figure 9.23, the 
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attenuation is obviously smaller than —11 dB. MATLAB program used 

w=2 0; 

while 20*logl0(100*abs(j*w+l0)/ 

abs( ( (j*w) ~ 2 + 2 * j *w+2)*(j*w+20) ) )<-11; 

w=w-l; 

end 

This program produces uj C gnew = 10 rad/s. In steps 5 and 6 the phase-lag-lead 
controller zeros and poles are obtained as —p c 2 = —17.3205, — 2+2 = —5.7735 for 
the phase-lead part and —p c \ = —0.1925, — z c \ = —0.5774 for the phase-lag part; 
hence the phase-lag-lead controller has the form 

s + 0.5774 w s + 17.3205 
~ s + 0.1925 X s + 5.7735 


Gm=lnf dB, (w= NaN) Pm=56.34 deg. (w=4.738) 



Frequency (rad/sec) 


Figure 9.24 Bode diagrams of the compensated system 
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It can be seen that the phase margin obtained of 56.34° meets the design requirement 
and that the actual gain crossover frequency, 4.738 rad/s, is considerably smaller 
than the one predicted. This contributes to the generally accepted inaccuracy of 
frequency methods for controller design based on Bode diagrams. 

The step responses of the original and compensated systems are compared in 
Figure 9.25. The transient response of the compensated system is improved since 
the maximum percent overshoot is considerably reduced. However, the system 
rise time is increased due to the fact that the system bandwidth is shortened 
(ucgnew = 4.738 rad/s < uj cg = 7.668 rad/s). 



Figure 9.25 Step responses of the original (a) and compensated (b) systems 


o 
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9.5 MATLAB Case Study 

Consider the problem of finding a controller for the ship positioning control system 
given in Problem 7.5. The goal is to increase stability phase margin above 30°. 
The problem matrices are given by 



—0.0546 

0 

0.5435" 


0 

A = 

1 

0 

0 

, B = 

0 


0 

0 

-1.55 


1.55 


C = [0 1 0], 


D = 0 


The transfer function of the ship positioning system is obtained by the MATLAB 
instruction [num, den] =ss2tf (A, B, C, D) and is given by 

= 0-8424 

v ’ s(s + 1.55)(s + 0.0546) 

The phase and gain stability margins of this system are Pm = —19.94° and 
Gm = —15.86 dB, with the crossover frequencies u cp = 0.2909 rad/s and uj cg = 
0.7025 rad/s (see the Bode diagrams in Figure 9.26). From known values for the 
phase and gain margins, we can conclude that this system has very poor stability 
properties. 

Since the phase margin is well below the desired one, we need a controller which 
will make up for almost a 50° increase in phase. In general, it is hard to stabilize 
systems that have large negative phase and gain stability margins. In the following 
we will design phase-lead, phase-lag, and phase-lag-lead controllers to solve this 
problem and compare the results obtained. 
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Gm=-15.86 dB, (w= 0.2909) Pm=-19.94 deg. (w=0.7025) 



Frequency (rad/sec) 


Figure 9.26: Bode diagrams of a ship positioning control system 


Phase-Lead Controller: By using Algorithm 9.1 with (f) max = 50° + 10° = 60° 
we get a phase margin of only 23.536°, which is not satisfactory. It is necessary 
to make up for <f> ma x = 50° + 27° = 87°. In the latter case the compensator has 
the transfer function 


Gc(s) = 76.31 


s + 0.2038 


s +15.55 
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Figure 9.27: Bode diagrams for a ship positioning system: 

(a) original system, (b) phase-lead compensated system 

The gain and phase stability margins of the compensated system are found from 
the above Bode diagrams as Gmc = 15.1603 dB, Pmc = 30.0538°, and the 
crossover frequencies are Lo cpc = 1.7618rad/s, uj cgc = 4.6419rad/s. The step 
response of the compensated system exhibits an overshoot of 45.47% (see Figure 
9.28). 


Phase-Lag-Lead Controller: By using Algorithm 9.3 we find the compensator 
transfer function as 


G c (s) = 7.524 " + °- 15 " x 0.1329 " + °- 016 

L \ / „ i i ono 


s + 1.203 


s +0.002125 
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Figure 9.28: Step response of the compensated system with a phase-lead controller 




Figure 9.29: Bode diagrams for a ship positioning control system: 
(a) original system, (b) phase-lag-lead compensated system 
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The phase and gain margins of the compensated system are given by 
Pmc = 39.6694°, Gmc = 14 dB and the crossover frequencies are uo C gc = 
0.4332 rad/s, u cpc = 1.2401 rad/s. 

From the step response of the compensated system (see Figure 9.30), we can 
observe that this compensated system has a smaller overshoot and a larger rise time 
than the system compensated only by the phase-lead controller. 



Figure 9.30: Step response of the compensated system with a phase-lag-lead controller 

Phase-Lag Controller: If we choose a new gain crossover frequency at u) cgnew = 
0.03 rad/s, the phase margin at that point will clearly be above 50°. Proceeding with 
a phase-lag compensator design, according to Algorithm 9.2, we get |AG(j'0.03)| = 
290.7390 and a = 0.034, which implies z c = 0.003 and p c = 1.0319 x 10 -5 . Using 
the corresponding phase-lag compensator produces very good stability margins 
for the compensated system, i.e. Gm = 32.91 dB and Pm = 54.33°. The 
maximum percent overshoot obtained is much better than with the previously used 
compensators and is equal to MPOS = 18%. However, the closed-loop step 
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response reveals that the obtained system is too sluggish since the response peak 
time is t p = 95.2381 s (note that in the previous two cases the peak time is only 
a few seconds). 

One may try to get better agreement by designing a phase-lag compensator, which 
will reduce the phase margin of the compensated system to just above 30°. In order 
to do this we write a MATLAB program, which searches the phase Bode diagram 
and finds the frequency corresponding to the prespecified value of the phase. That 
frequency is used as a new gain crossover frequency. Let Pm = 35° = 0.6109 rad. 
The MATLAB program is 

w=0.1; 

while pi+ 

angle(1/((j*w)*(j*w+l.55)*(j*w+0.0546)))<0.6109; 

w=w-0.01; 

end 

dG=0.8424 *abs(1/((j*w)*(j*w+1.55)*(j*w+0.0546) )); 

This program produces uj cg new — 0.07rad/s and |AG(j'0.07)| = 87.3677. From 
step 4 of Algorithm 9.2 we obtain the phase-lag controller of the form 



s + 0.009 
s + 0.000081 


x 0.0114 


The Bode diagrams of the compensated system are given in Figure 9.31. 
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Gm=23.48 dB, (w= 0.2651) Pm=31.44 deg. (w=0.06043) 



Frequency (rad/sec) 


Figure 9.31: Bode diagram of the phase-lag compensated system 


It can be seen that the phase and gain margins are satisfactory and given by 
Pm = 31.44° and Gm = 23.48 dB. The actual gain crossover frequencies are 
u'cynttr = 0.06043 rad/s and u cpnBW = 0.2651 rad/s. 

The closed-loop step response of the phase-lag compensated system, given in 
Figure 9.32, shows that the peak time is reduced to t p = 50.15 s—which is still fairly 
big—and that the maximum percent overshoot is increased to MPOS = 45.82%, 
which is comparable to the phase-lead and phase-lag-lead compensation. 
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Figure 9.32: Step response of the phase-lag compensated system 

Comparing all three controllers and their performances, we can conclude that, 
for this particular problem, the phase-lag compensation produces the worst result, 
and therefore either the phase-lead or phase-lag-lead controller should be used. 


o 




